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Abstract 

A system of resonance atoms is placed in a medium with developed 
polariton effect. In the spectrum of polariton states there can exist a 



Ch , band gap. If an atom with a resonance frequency inside the polariton 



gap is incorporated into the medium, the atomic spontaneous emission can 
be strongly suppressed. This is what is called the localization of light. 



H ' Nevertheless, an ensemble of resonance atoms inside the polariton gap can 

radiate if their coherent interaction is sufficiently strong. Conditions when 
coherent polariton radiation can appear and the properties of this coherent 
polariton emission are studied. 



1 Light Localization 

The phenomenon of hght locahzation appears in three-dimensional periodic di- 
electric structures, in which, due to periodicity, an electromagnetic band gap 
develops. Then spontaneous emission with a frequency inside the band gap can 
be rigorously forbidden [1-3], because of a severe depression in the photon den- 
sity of states for those frequencies which remain in the spectral gap between the 
upper and lower branches. This kind of samples, in which photon band gap de- 
velops because of the structure periodicity in real space, has been called photonic 
band-gap materials. The appearance of gaps in the spectrum of photon states, 
due to real-space periodicity, is similar to the formation of gaps in the spectrum 
of electrons in a periodic lattice potential [4]. 

If resonance atoms are incorporated in a band-gap material, so that their 
transition frequency is inside the gap, then the effect of light localization [2,3] 
can arize. To formulate explicitly what this effect means, let us consider the 
average 

s{t) =< a'{t) > 

of the population-difference operator a^. The average here implies the statistical 
or, at zero temperature, quantum-mechanical average. Under the localization of 
light one understands [5,6] that 

lim s{t) > -1 . 

The light localization becomes possible due to the formation of photon-atom 
bound states [5-7]. If a collective of identical impurity atoms is incorporated 
into a medium with a photon band gap, so that their transition frequency is 
inside the gap, and their spacing is much less than the transition wavelength, 
then a photonic impurity band is formed within the photonic band gap [7]. Elec- 
tromagnetic coupling of neighboring atoms takes place by means of an effective 
resonance dipole-dipole interaction. If this interaction is sufficiently strong, then 
electromagnetic radiation can propagate inside the impurity band [7]. 



The formation of photon band gaps in photonic band-gap materials is similar 
to the well known polariton effect of the appearance of photon bands due to the 
interaction of light with collective excitations of dense medium [8,9]. Physical 
processes are, actually, the same in both types of materials. The difference is 
only in the nature of scatterers which light interacts with. In artificial photonic 
band-gap materials, a suppression of the photon density of states over a narrow 
frequency range results from multiple photon scattering by spatially correlated 
scatterers. In natural dense media, such as dielectrics or semiconductors, a fre- 
quency gap for propagating electromagnetic modes develops as a result of the 
photon interaction with optical collective excitations, such as optical phonons, 
magnons, excitons, and so on. Photons in a medium, coupled with collective 
excitations, are called polaritons. 

When a single resonance atom is placed in a frequency dispersive medium 
whose polariton spectrum has a gap, and the atomic transition frequency lies 
inside this gap, then a polariton-atom bound state appears causing a significant 
suppression of spontaneous emission [10,11]. The physical picture explaining 
this suppression is as follows. Let us imagine an atom in a medium, with the 
atomic transition frequency within the polariton gap. If this atom is initially 
excited, then it tends to become deexcited emitting a photon. However, since 
the propagation of photons inside the polariton gap is prohibited, the emitted 
photon is scattered, by collective excitations, back and is again absorbed by the 
atom. Thus, the atom cannot get rid of a photon and is doomed to stay excited. 
Conversely, if the atom is initially in the ground state, it continues to be in 
that state, since there are no photons around to excite it. The supression of 
spontaneous emission of an atom is termed localization of light. As is explained 
above, the effect of light localization can be expressed as the inequality s(t) > —1 
for the average population difference, valid for all times. 

The formation of polariton-atom bound states has been studies for a sta- 
tionary case [10,11]. In dynamical picture, the population difference of an atom 



satisfies the equation 

ds 

in which 71 is a level width and (, a stationary value of the population difference 
defined by a solution to the stationary problem. From the dynamical equation 
one has 

s{t) = (so - C)e-''' + C , 

where Sq = s(0). If the stationary value ( > —1, then Mmt-^^ s{t) = ( > —1, 
which implies the localization of light. The complete suppression of emission 
corresponds to C = sq] then s(t) = (. Note that the linewidth 71 is caused by 
vacuum quantum fiuctuations and is always nonzero, irrespectively what medium 
the atom is placed into. 

If a collection of resonance impurity atoms is doped into a medium with a 
polariton band gap, then, in the same way as for photonic band-gap materials 
[7], an impurity band can be formed within the polariton gap [12,13]. Then 
electromagnetic radiation can propagate in such an impurity band. In order that 
such an impurity band be formed, the spacing of resonance impurity atoms in 
the medium should be much smaller than the radiation wavelength. If it is so, 
then for a group of atoms a sufficiently strong effective interaction, caused by 
photon exchange, can develop. This interaction collectivizes the atoms that start 
radiating coherently [14]. In this way, the suppression of spontaneous emission 
for a single atom can be overcome by a group of atoms radiating coherently. 

The situation when a single atom cannot radiate inside the polariton gap 
but a collective of strongly interacting atoms can radiate reminds the following 
related case. If a sample with a polariton band gap is irradiated by a monochro- 
matic electromagnetic wave with a frequency within the polariton gap, then the 
incident light cannot propagate through this medium because of total reflection. 
However, if the incident intensity is large enough, the light can penetrate into 
the dense media even when propagating inside the polariton gap [15,16]. In such 



a case, analogously to that of coherent radiation, the possibility of the radiation 
propagation inside the band gap is due to nonlinear effects. 

Here and in what follows we use the term "atom" in the general sense, 
implying under a resonance atom any two-level object. Depending on radia- 
tion frequencies, this could be atoms as such, molecules, nuclei, or quantum 
dots and wells. The latter case is of special importance for semiconductors. 
Really, the polariton effect is well developed in many semiconductors, for in- 
stance, in CuCl, CuBr, CdSe, ZnSe, GaAs, GaSb, InAs, AlAs, SiC. 
The characteristic frequencies, where the polariton band gap arises, are as fol- 
lows [4] (see also [12-14]). For example, in GaAs the polariton gap having 
the width A = ri2 - ^i = 4 x W^s'^ lies between Vti = 5.1 x lO^^s"^ and 
^2 = 5.5 X lO^^s"^; the linewidth being 7i/r2i = 1.2 x 10"^. In SiC the polariton 
gap A = 3 X lO^^s"^ is between ^i = 1.5 x lO^S"^ and 1^2 = 1-8 x lO^^s"^; with 
the linewidth 7i/r2i = 3 x 10^^. In all cases, for the relaxation parameter 72 one 
has 72/^1 ~ 10~^. As is seen, the polariton band gap in such semiconductors is 
located in the infrared region. Therefore, resonance radiation for this region of 
frequencies could be presented by quantum dots and wells. Keeping in mind the 
feasibility of different radiating objects, we continue, for the sake of simplicity, to 
use the term "resonance atoms". 

2 Basic Equations 

The total Hamiltonian is given by the sum 

H = Ha + Hf + H„i + Haf + Hmf , (1) 

consisting of atomic. Ha-, field, Hf, matter, Hm-, atom-field, Haf-, and matter- 
field, Hmfi Hamiltonians. In the atomic Hamiltonian 

1 " 
Ha = -Y.^,{l + al) (2) 



the index i enumerates the atoms, ujq is a transition frequency, and af is a pop- 
ulation difference operator. Here and in what follows we set h = 1. The field 
Hamiltonian 



^ Stt 



,2 ,2 

E (r)+ H (r] 



d7 (3) 



has the standard form in which E is electric field and H=V x A is magnetic field, 
with a vector potential A satisfying the Coulomb gauge condition V ■ A= 0. The 
Hamiltonian of matter represents optic collective excitations and can be modelled 
by an ensemble of oscillators, 

N' -^"^ 1 Af' 3 



^"^ = ES + 2EE^r<<' (4) 

i=l ij af3 



where the index i = 1,2, . . . ,N' enumerate lattice sites. Pi and Mj are momentum 



and displacement operators, respectively, and D""- is a dynamical matrix. The 



a/3 



atom-field interaction is described by the Hamiltonian 

1 ^ -> - ^ 
Haf = --T.Ur^)A{r.), (5) 

'^ 1=1 
which corresponds to the dipole approximation with the transition current 



Ja (ri) = iujQ I (7+ d -(^i d] , (6) 



where af and a^^ are the rising and lowering operators, respectively, and rf is a 
transition dipole. The matter-field interaction can be presented as 

1 AT' ^ _^ 

^m/ = — E7m(r,)A(r,), (7) 



with the matter current 



^ (^ \ — ^ ^ 

dm {' i ) -r J , 

m 



in which e and m are charge and mass, respectively. 

The commutation relations for the operators introduced above are 

[at, a-] = 5^jat , [al af] = ±25,j(jf , 



[E»(7), A" {7 )] = t47rc6apS{7 -7). 

Using these relations and the Heisenberg equations of motion, we get the Maxwell 
operator equations 

1 9 i - - 47r - 



Id A 
c dt 


= -E 


with the total density of 


current 



where the notation 



N N' 

Jir)=J2 Ja i7,)5i7 - r,) + ^ J™(r,)5(? - ?,) . (10) 



-iujoo-i + koa- d ■ Ai , (11) 



^^ 2A;o((T+ d +^r d)- Ai , (12) 



Ai=A(7i,t), ko = — 

c 



For the atomic variables 


>, we 
da- 


find 




dt 




and 


daf 





is used. From (9), with the Coulomb gauge condition, we have the wave equation 

whose solution reads 

J (?, t) =i (r , t) + - I nr ,t-Jr-r \/c) ^ ^, ^ ^^^^ 

where Av, being a solution of the related homogeneous equation, corresponds to 
vacuum fluctuations. Substituting the density of current (10) into (14) yields for 
the vector potential at the point r j the expression 

Ai (t) =Av (7i, t)+ Aa {7i, t)+ Am (7i, t) , (15) 

in which the first term is caused by vacuum fiuctuations, the second term, 

(16) 



Aa {7i, t) = iko J2 



V(*-^)5'-7(*-^U 



j(^i) ^ij 



where 



T- ■ ^ T ■ ■ 



/y ^^ /y /y 

' ij — ' i j 1 



is a vector potential generated by radiating atoms, and the last term, 



1 ^ 1 



Amir,t) = - J2 — Jrn (rj,t- ^ 



C j(^i) ^iJ 



(17) 



is due to local electric currents in the medium. In the vector potentials (16) and 
(17) the self-action parts are excluded. Instead, we shall add to Eqs. (11) and 
(12) the terms describing the level width and the line width. 



2 1 

71 = ^hdo = Y ' 



72 



1 



where do = \ d \- l^ this way, introducing the effective electric induction 



D^ (t) = ko 



Av (r,t)+ Am (r,t) 



we come to the equations 

dar 



dt 



[iuJo + 72)(Ti +<yi d ■ Di + 



N 



+ ikl<Jl d- E — 



ji^i) ^iJ 



^ij 



''[t-f)" -7(,*-f)'i 



Ti 



(19) 



and 



dal 
dt 



-7i(^'-C) + 2(af rf +a- d)- Di 



N 



i2kl{at d +a- d)Y.— 



ji^i) ^iJ 



^tit-"-^]I-a7(t-^-^]d 



(20) 
c / •' \ c / ] 

The retardation effects in these equations can be treated in the quasirela- 
tivistic approximation. This means the following. In the nonrelativistic limit, 
when c ^ oo and /cq -^ 0, from (19) would follow a^ ~ exp{—iujQt). In the 
quasirelativistic approximation, we set 



ar(t 



o-j (t)exp{ikori 



(21) 



Define the statistical averages 



Ui =<(Ti > , Si =< al > (22) 



over atomic degrees of freedom. Tlien from (19) and (20), in the semiclassical 
approximation, we obtain 



dt 
and 



dzL ' — ^* — ^ Q — * V— -> — ^* — * 

-^ = -{iuo + 72)Mi + Si{di ■ Di) + ikpi d ■ Yl i'^ij'^*j d -^ijUj d) (23) 



where 



ds' — ^* — ^ — ^ 

-^ = -lii-Si - C) - 2« d +u, d)- A - 

t2kliu* d +u,d)-J2 (^*i«j d -VijU, d) , (24) 



^ exp(?A;orij) 

'^ij — 



The semiclassical approximation is a kind of the mean-field approximation. 
In the spirit of these, we may make the following mean-field approximation 

N 
Y^ ipijUj ^ UiY Vij = Ui(pi , 

where 

J^ J^ exp(?A;orij) 



The factors (pij and (pi describe local fields. 
Introduce the local-field shift 

A, = ^,g's , g'^3^J2 ^^T^^ ' (25) 

also called the cooperative Lamb shift [17], and the effective atom-atom coupling 

parameter 

^ ^44 ^ sin(fcor^j) 

9 — ^, 2^ u ^ ■ l^"J 

72 j^^i) Kor-ij 



These quantities enter into the definitions of the effective radiation frequency and 
radiation width, 

n = uo + Ai^, T = 72(1 - gs) , (27) 

respectively. 

Involving these notations and keeping in mind that 

Ui = u(ri,t), Si = s(ri,t), Di=D (ri,t), (pi = (p{ri), 



we transform Eqs. (23) and (24) to the form 
du 



dt 



-{in + T)u + s d ■ D+ik^sip*u*{d 



(2J 



and 



— = -472^|mP-7i(s-C)-2(m* rf +ud)-D 



tZithf] 



(29) 



^*{u* d f-v{udf 

Since m is a complex variable, we have to supplement Eqs. (28) and (29) by an 
equation for either u* or \u\^. For instance, for \u\^ we get 

d\u\'^ 



dt 



-2V\u\^ + s{u' d +u d)- D + 



+ ikls i^*{u* d f -^{u df 
Equations (29) and (30) give for the Bloch vector the equation 
d 



(30) 



dt 



(s^ + 4|m|2) = -872|Mp - 27i(s - C)s 



The derived equations (28), (29), and (30) are the basic equations describing 
nonequilibrium processes in the system of resonance atoms interacting with po- 
lariton field. 



10 



3 Scale Separation 

Equations (28), (29), and (30) can be solved using the scale separation approach 
[18-20]. To start with, we need to define what small parameters we have. 

The standard small parameters are related to the relaxation parameters 71 
and 72, for which 






^«1. 



It is reasonable to suppose that 



n 



< 1 



<i, 



(31) 



(32) 



although r can become much larger than 72. Assume also that 



d ■ D 



n 



<i, 



(33) 



which means that the interaction of atoms with matter does not change drastically 
the properties of the atoms. Under the validity of small parameters (31) to (33), 
the variable u has to be considered as fast, compared to s and j-up that are to be 
treated as slow. Accepting the variables s and |-up as quasi-integrals of motion, 
we keep them fixed when solving Eq. (28). Then the solution for the fast variable 
is 

u{t) =uoGi{t) +u*QG2{t)+ 

+ sl* f [Gi{t - r) + G2{t - t)] D {T)dT , (34) 

JO 

where uq = m(0) and the Green functions are 

Ai — a* \ N . I X2 — a* 



G,{t) 



.Mt 



Ai — A2 



Ai — A2 



,Ait 



^\2t 



K 



G2{t) = 

Ai — A2 

-(zf] + r), h = -ikls^*{df 



11 



Taking into account the existence of small parameters, we may write 

and (34) can be simplified to 

u{t) = e-(^^+r)* L + sl [' e(^^+r)- D {r)di 

Jo 



(35) 



The found fast variable (35) is to be substituted into the equations (29) and 
(30) for the slow variables and the right-hand sides of these equations are to be 
averaged over time and over the degrees of freedom corresponding to collective 
excitations of matter [18-20]. Recall that the quantities in (22) were defined as 
the averages over atomic degrees of freedom. The double averaging, over time 
and over the matter degrees of freedom, for a function F{t), depending on time 
and on the operators of collective excitations, is defined as 

« F »= lim - r < F(t) > dt , (36) 

r^oo T Jo 

where the angle brackets imply the statistical averaging over the matter degrees 
of freedom. The usage of the same angle brackets for denoting the statistical 
averaging over the atomic and over matter degrees of freedom should not yield 
confusion, since at the present stage the atomic degrees of freedom do not arise 
being averaged out earlier. Therefore, in the definition (36) and in what follows 
the statistical averaging always concerns only the matter degrees of freedom. 
Let us introduce the parameter 



a =<< 



^^Tt f\i,n+T)r 2*. ^ {^r)dr 
Jo 



» (37) 



characterizing the strength of interaction between the atoms and matter. Thus, 
the quantity (37) can be called the atom-medium coupling parameter. 

When substituting the fast variable (35) into the equations (29) and (30) for 
the slow variables and averaging, according to (36), the right-hand sides of the 



12 



latter equations, we may notice the following useful property. The fast variable 
(35) can be written as the sum u = Ui + U2 oi the term 

Ml = Mo exp { — {iQ + r)t} , 

not depending on the field D of matter, and of the term 

-(iQ+r)< ^* /■* (if7+r)r 



u, = e-^^'^+'^'s d / e^^^'+^> D {T)dT , 
Jo 

depending on the matter field. Define the function 

w = \ui\ = \uo\ e 

By this definition, the function w must satisfy an equation that follows from the 
equation for |-up where the matter field D is set zero. For |-up we have 

|-Up = W + ulU2 + U^Ui + |-U2p • 

When averaging, according to (36), we take into account that 

<< ulu2 + U^Ui »= . 

This is because of two reasons. First, the terms ui and U2 oscillate, in general, 
with different frequencies. Second, the term U2 is a linear combination of op- 
erators of collective excitations in matter. In this way, averaging |-up over fast 
variables, we get 

|Mp = W+ « |M2p >> , << |M2p >>= tt-S^ • 

This consideration suggests that it is convenient to introduce the slow variable 

w = \u\'^ — as'^ , (38) 

for which the evolution equation should have a form simpler than for |-up. Really, 
averaging the equations (29) and (30) for the slow variables, we obtain 

ds 

^ = -%72W-7i(s-C) , (39) 

13 



where the transformation (38) is used, and 

^ = -2^^{l-gs)w. (40) 

From these two equations one can derive one equation 

d s ds 

— + (2 + 7 - 2gs)- - 2^gs' + 27(1 + gC)s - 27C = , 

in which 7 = 71/72 and time is measured in units of 7^"*^. 



4 Models of Matter 

Before analysing equations (39) and (40), let us consider some examples defining 
concretely the matter field D- Suppose, first, that the matter consists of a set of 
random scatterers, such that 

rf ■ B= e , (41) 

where ^ is a stochastic field defined by the averages 

< e >= , < ler >= 7' • (42) 

Then the coupling parameter (37) is 

As another example, consider the matter modelled by the white noise, when 

1 ■ D= m , (44) 

where the white-noise stochastic variable C,{t) is defined by the averages 

< m >= , < c{t)m >= 275(t - 1') , (45) 

where the angle brackets mean a stochastic averaging. Then for the coupling 
parameter, we get 

« = ^ • (46) 

14 



In the third example, we model the matter by an oscillator, so that 

1-D=^ {h^e~'^' + hie'^') , (47) 

where 6^^ and h]^ are the annihilation and creation operators satisfying the Bose 
statistics, and for which the statictical averaging gives 

< h^hl, >= 1 + n^ , 



< 6^,6^; >= n^ , 



(4J 



with riu) being an effective weight of excitations of a frequency uo. Then the 
coupling parameter (37) is 

_ I |2 '']2£, I '- '^ f^u) (AQ\ 

"~'^' [isi-ujy + v' iyt + uy + T^l' ^ ' 

The strongest coupling between the impurity atoms and matter happens at the 
resonance, when oj = Vt, and a = n^|7/rp. 

Finally, we consider a more realistic situation when the effective electric in- 
duction of matter is defined by the relations (18), (17), and (8), so that 



D^ 

with the momentum operator 

P^ it) 



eko X 1 



E 



Pi t 



mc .^_^.^ Tij 



u 



(50) 






ks 



V 2N' 



6fcsexp{?(A; ■ Tj -Ukst)} - bl^exp{-i{k ■ rj -Ukst)} 



(51) 



in which Uks = ^-ks is a spectrum of collective excitations; k being a wave vector; 
s = 1,2,3, a polarization index; e^^ is a polarization vector; A^' is the number 
of lattice sites. The annihilation and creation operators of collective excitations 
satisfy the Bose statistics and have the following statistical averages 



< bj^gbk's' >— nksSkk'Sss' , 



< bksbk's' >— . 



In this case, for the coupling parameter (37), we obtain 

nks , 1 + nks 



_ kor-e >-^ 

tt — -TTJJ 2^ Jkslks^ks 



2N' 



ks 



{n-ujksy + T^ {n + ujksy + T\ 



(52) 



(53) 
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where 



and 



Iks = kl\ d ■ 6/ 



Jks — 



r^. = 



rac^ 



_'£ks 



N' exp ji ( /c ■ r^ , ^ , .^ 



(54) 



(55) 



Again, it is clear that the couphng parameter (53) is the most strongly influenced 
by resonance collective excitations with uj^s ~ fi. 



5 Transient Regime 

Consider the times shorter than 7f^. Then the term containing 71 in (39) can be 
omitted. In this case, using the second relation from (27), we have 



dr 



Ag^-f^w 



dw 



dt dt 

These two equations can be reduced to one, 

d^r dr 



-2rw 



(56) 



integrating which we get 



dT 



+ r' = 7o, 



7o being an integration constant. The last equation is a Riccati equation whose 
solution is 



r = 7otanh ( , , 



7o = — , 

To 



where to is another integration constant. Using relation (27) gives 



7o , . ft-to 
tanh 



9l2 V T-Q 

The first equation in (56), together with (57), yields 

't-to 



(57) 



(5J 



w = . ° „ sech^ 



2^,2' 



%'72 



To 



(59) 



16 



And from (38) we find 



\u\ 






2o,2'- 



Ag^Y2 



as 



(60) 



The integration constants 70 and to are to be defined from the initial conditions 



m(0) = Mo 



s(0) = So • 



(61) 



From the latter we obtain for the effective radiation width 70 the expression 

2 1^2 I /I 2 2 /i |2 2\ 

7o = ro + 4^ 72 (foI -ttoSo > 



(62) 



in which a^ is a ai t = 0, when s = sq, 



To = 72(1-5-50) 



and the delay time 



To 



In 



7o 



7o + ro 



The radiation width can be written in the form 

70 = 2|5(So|72(«c-tto)^^^ , 

where the critical atom-matter coupling parameter 

(l-5So)2 + 4(7>oP 



ar 



2„2 



Ag^s'o 



(63) 



(64) 



(65) 



is introduced. Expression (64) is a direct consequence of (62) for all g and a^. 

It is necessary to stress that the atom-matter coupling parameter (37) cannot 
surpass the critical value (65). If this would happen, then the radiation width 
(64) would become imaginary and, instead of (57), we would have 

r = -|7o|tan - — - , |t-o| = -^ — r , 

V rol / |7o| 

l7o| = 272|fi'So|Vao - "c (ao > "c) • 

The delay time (63) would be 



to = |To|arctan 



l-fi-so 



2|5fSo|V"o - ttc 
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And for the solutions (58) and (59) we would get 
s = SoSgn{gsQ)y/ao — Cretan 



t-to\ 1 



To y 9 

2( \ 2 (^ ~ ^^ 

w = — So(ao — «c)sec 

V To 

The effective width F, as well as the solutions s and w, become divergent at 

tn = to + ^(l + 2n)|ro| (n = 0, 1, 2 . . .) • 

Certainly, this behaviour is unphysical and it means that some conditions, under 
which the method of scale separation has been used, are, probably, not valid any 
more. This is really the case since, when F and s diverge, conditions (32) do not 
hold true. Consequently, when a^ > ac, we cannot separate solutions onto fast 
and slow, all of them oscillating equally fast. At the same time, the existence 
of slow solutions is a characteristic feature of developed coherence. Thus, the 
absence of slow solutions suggests that coherence cannot emerge in the system. 
From the physical point of view, all this sounds quite understandable. There 
should be a threshold for the strength of interactions of atoms with matter af- 
ter which such strong interactions destroy the correlation between atoms, thus, 
destroying their coherence. In this way, the inequality 

ao < etc (66) 

is a necessary condition for the applicability of the scale separation approach and, 
at the same time, a condition for the possibility of coherent radiation of doped 
atoms. 

The maximal level of coherence develops at the time t = to when 

s(to) = - , w(to) = sl{ac - ao) , \u(to)\^ = {ac - ao)sl + — . (67) 

9 9 

For the times much longer than to, Eqs. (58) to (60) give 

s ^ - f 1 - -) (t » to) , 

9 V 72/ 
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«;~4Sexp(-27ot) , (68) 

9 l2 9 \ 12/ 

However, the asymptotic behaviour given by (68) is vahd only for t <^Ti. 

Consider the case when both the atom-matter and atom-atom couphng pa- 
rameters are small, i.e. 

ao ^ ttc , Ifi*! ^ 1 • (69) 

Using the first of inequalities in(69), we have from (64) 

70 - 72 [(1 - 9-%? + %>o|'] '^' (l - ^) • 
The latter expression, with the second inequality in (69), reduces to 

7o -72(l-^So) (l --^j . 
The critical parameter (65) becomes 

a,~(Vs2)-i {g<^l). (70) 

Employing this, we find 

70 ^ 72(1 - 9Sq - ^aog^sl) , (71) 

valid for small coupling parameters as in (69). For the delay time (63), we get 

to^^^\nao9'sl , (72) 

272 

which tends to —00 if either ao or 9 tends to zero. This implies that, under 

conditions (69), an essential level of coherence does not evolve. 

Let us analyse the case when 

ao ^ «c , Ifl*! > 1 • (73) 

Then the critical parameter (65) is 

a^^^(sl + Muo\'-2'-^] . (74) 

4so \ 9 J 
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The radiation width (64) becomes 

i2\g\ 



7o 



4|mo| 



9 A \ 1 9 r\ 9 

So + 4|mo| -2aoSo 



(75) 



with the corresponding radiation time 



To 



To 



g\Js^ + 4\uo\ 



(76) 



For the delay time (63), we find 



T"0 



\g\{sl + 4|moP - 2aoSo) + 9SoJ4 + M^o 



(77) 



\g\{sl + 4|moP - 2aoSo) - 9So\Jso + 4|mo|2 

If the process develops from an initially incoherent state, when uq = 0, then the 
radiation width (75) is 



7o - 72|fi'So| 1 - 2ao (^o = 0) . 

V 9-30 J 



(7J 



Thence, the delay time (77) becomes 



— ~\ r ^^ 

2|fi'So| 



1 — 2ao + e 



1 — 2an — e 



(79) 



where 



e = sgn{gso) = ±1 . 
As far as for Uq = and \g\ ^ 1, the critical parameter (74) is 

"c^ ^ (1^1 > 1, Mo = 0) , 
then the inequality ao <^ ac implies ao ^ 1. Hence, we may simplify (79) as 

T2 



tn 



25-50 



llnctol . 



^0) 



After the time (80), the population difference tends to 

1 + e 



-eso(l - 2ao) + 



9 



it > to) 
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li g > 0, then esq = \sq\, while for g < 0, one has ssq = — |so|- Combining both 
these cases, we get esq = sgn{g)\sQ\. 

Assume now that the atom-matter couphng parameter is close to its critical 
value (65), but the atom-atom coupling is arbitrary, 

^^^^^ « 1 (V^) . (82) 

Oic 

Then the radiation width (64) tends to zero, as ao ~^ o^c, and respectively, the 
radiation time Tq = 'Jq^ tends to infinity. For the delay time (63), we find 

to^jYZ N^("c-ao)^ , (83) 

while the radiation time is 



T, 



2 



^\gso\ 

When ao -^ a^ then to -^ 0, and for the functions (58) and (59) we have 



s ~ 
9 



- - 2|so|sgn((7)Vac-ao (l - 26-^^°*) , 

w ~ Asl{ac - ao)e"^^"* (t > to) . (85) 



There is a suppression of self-organized coherence in the system of atoms, their 
radiation being almost completely due to the pumping by matter excitations 

1 n I 12 '^ 

S ^ - , W ^ , \u\ ~ ^7 • 

9 9 

Although the coherent relaxation may happen provided that To ^ T2, that is, 

\gsQ\{ac-aof^ > 1 , 

which corresponds to superradiant emission. 

Note that analysing the properties of the solutions to equations (39) and 

(40), we talk about radiation processes keeping in mind the following. The total 

radiation intensity of atoms can be approximately defined in the usual way as 

ds 
lit) = -Nfiu^— . (86) 
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For a more accurate definition of radiation intensity see e.g. Refs. [21,22]. From 
(86), using equation (39), we find 

J(t) = hoh{t) + hncit) , (87) 

where the first term 

hoh{t) = ANghuJo^2W (88) 

had the meaning of the coherent radiation intensity, and the second, 

I^ncit) = NhuJo-fi{s - , (89) 

corresponds to the intensity of incoherent radiation. The latter is always propor- 
tional to the number of atoms A^, while the radiation intensity (88) is proportional 
to Ng. For a concentrated sample, whose linear size is much smaller that the ra- 
diation wavelength, we have g ^ N, and the radiation intensity (88) becomes 
proportional to N^, which is typical of superradiance. In this way, the solutions 
s and w define the temporal behaviour of the incoherent radiation intensity (89) 
and of the coherent radiation intensity (88), respectively. For instance, using the 
solution w given by (59), with the radiation width (64), we obtain the intensity 
of coherent radiation 

2/ N 1,2 /^ ~ ^0 



Icohit) = 4NghuJoi2st{ac - ao)sech' ( — ^ ) • (90) 

The latter shows that, if ao — ^ «c, then Icoh -^ 0. 

6 Close— to— Stationary Regime 

In the previous section the transitient regime is considered corresponding to times 
t <^Ti. For the times comparable or larger that Ti, we cannot neglect any more 
the term with 71 in equation (39). In the intermediate stage, when t ~ Ti, an 
exact solution of Eq.(39) and (40) is not available. Here we have to resort to 
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numerical calculations, which will be the subject of a separate paper. But it is 
possible to give an analysis for asymptotically large times, when t ^ Ti. 

The following analysis assumes that g ^ 0. Since, ii g = 0, the solutions to 
Eqs.(39) and (40) are 

s = C + {so- C)e-''' , 

w=(|Mo|'-«So)e-'^'* (^ = 0), (91) 



which describes the relaxation process of a single atom. In such a case, if there 
is the localization of light, then (^ = sq, and (91) gives s = sq. 

If A^ impurity atoms are doped into the matter, then g ^ 0. The resonance 
dipole-dipole interactions of a pair of atoms with a transition frequency inside 
the photon gap have been studied in several works [7,23,24]. The conclusion of 
these studies is that two closely spaced atoms, with transition frequencies in the 
gap, interact with each other by means of the virtual photon exchanges much 
in the same way as the atoms in vacuum. That is, if the atoms are separated 
from each other by a spacing much larger than the radiation wavelength, then 
each of them can be considered as a single atom. If the transition frequency of 
such an atom is inside the gap, then the phenomenon of light localization occurs. 
However, if the atoms are close to each other, with a spacing mush smaller than 
the radiation wavelength, then they practically do not experience the existence 
of the gap [7,23,24]. 

Consider the close-to-stationary regime, when t ^ Ti. Equations (39) and 
(40) can be written as 

with the right-hand sides 

Vi = -Ag-f2'w - -fi{s - C) , 

V2 = -2-f2{l - gs)w . (93) 
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Stationary points, or fixed points for Eqs.(92), (93), are given by tlie condition 
\/^ = V2 = 0. Tliis yields two stationary points: 



s* = C, wl = 



and 



1 



Wr, 



iiii-gC) 



(94) 



(95) 



' g ' 4725-2 

Tlie stability of these fixed points can be defined by the Lyapunov analysis. To 
this end, we need to find the eigenvalues of the Jacobian matrix 



J 



dVi dVi 
ds dw 



dV2 dV2 
ds dw 



(96) 



These eigenvalues are given by the expression 
1 



A^ 



71 + 272(1 - gs) ± (71 - 272(1 - gs)y - 32^^g'w 



1/2 



(97) 



Substituting here the values corresponding to the fixed points yields the Lyapunov 
exponents. For the stationary point (94), we have 



A| = -7i, Ar = -272(1 -50, 



(9J 



and for the stationary point (95), we find 



i± 



7i 



A? = -^ n ± 



1 + 8^(1-50 
7i 



nl/2' 



(99) 



The analysis of the Lyapunov exponents (98) and (99) shows that if 



K<1 



(100) 



then the fixed point (94) is a stable node, and the fixed point (95) is a saddle 
point. When 

K=l, (101) 
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both fixed points merge together becoming neutral, since Aj~ = A^ = 0. In this 
case, the system of equations (39) and (40) is structurally unstable. Equality 
(101) defines a bifurcation point. For the interval 

K ^C < 1 + ^ , (102) 

the fixed point (94) is a saddle point, while that (95) is a stable node. For all 
g( > 1, the point (94) is a saddle point. If 

^C > 1 + ^ , (103) 

072 

the stationary point (95) becomes a stable focus, since the Lyapunov exponents 
(99) take the form 

Af = -yT^cc;oo, (104) 



where 



_ 71 



80, 1 ^/2 

^(K - 1) - 1 

7i 



Suppose that for a single atom there occurs the localization of light, so that 
( = Sq. If many resonant atoms are doped into matter, but their interactions 
through the polariton exchange are not strong enough, so that qsq < 1, then 
the light localization prevails. This means that the light remains confined in the 
vicinity of the atoms. The confinement of light is demonstrated by the fact that 
the stationary point (94) is a stable node with si = Sq. But if the resonant 
interaction between the atoms is sufficiently strong, so that 

gso>l, (105) 

the deconfinement of light happens. Then the fixed point (95) becomes stable, 
while the point (94) looses its stability. The deconfinement of light is not com- 
plete, since s^ = 1/g < Sq, but it is partial. The portion of light that remains 
confined decreases with increasing g. The qualitative change of the asymptotic 
behaviour of solutions to a system of differential equations is called, in dynamical 
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theory, the dynamical phase transition. In our case, this happens ii g( = 1. The 
equahty gs^ = 1 separates the regions where hght is locahzed {qsq < 1) and where 
it is deconfined {qsq > 1). Therefore, the dynamical phase transition occurring 
at qsq = 1 corresponds to a transition that may be called the deconfinement of 
light or photon deconfinement. 

When the resonant interaction between atoms is so strong that inequality 
(103) holds true, then the stable stationary point (95) is a focus. This means 
that the solutions to the equations (39) and (40) display an oscillatory regime of 
motion when approaching the stationary point (95). Such an oscillatory motion 
is similar to that found for a concentrated system with the resonant frequency 
near the edge of a photonic band gap [5] and to that for two atoms with transition 
frequencies inside or shghtly outside a photonic band gap [24] . 

7 Coupling Parameters 

There are several characteristic quantities defining the behaviour of the system. 
These are the initial conditions Sq = s(0) and uq = u{0) and the coupling param- 
eters g, g', and a. Below we study the typical values of the latter. 

Recall that the coupling parameters g', defined in (25), and g, given in (26), 
have appeared in the evolution equations when treating the retardation effects in 
the quasirelativistic approximation (21). Without the latter approximation, we 
should deal with the integral-type equations [25] . Thus, the atom-atom coupling 
parameters g' and g describe the retardation or local-field effects. The values 
of these parameters essentially depend on the shape of the sample and on the 
spacing between atoms [17]. 

Accepting the equality 71 = 2kQdl/3, we have from (25) and (26) 

_ 371 y> sm{korij) ' _ ^ y^ cos^kprij) 

272 ^.(_,.) kor.j ' ^ 272 ^,(^^) kon^ 

If the radiation wavelength A = 27r/A;o is much smaller than the mean spacing, 
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a, between the atoms, then the sums in g and g' can be of any sign but with 
absolute values less than unity. As far as usually 71 <^ 72, the absolute values 
\g\ and l^f'l are small. If \g\ ^ 1 and \g'\ <^ 1, the impurity atoms almost do not 
interact with each other and their behaviour is practically the same as that of a 
collection of single atoms. 

In the opposite case, when A ^ a, the sums in g and g' can be estimated with 
a good approximation [17] as 

y. sm{korij) ^ it cosjkorij) ^ 3 

where p = N/V is the density of the doped atoms and it is assumed that A is less 
than the linear sizes of the sample in all directions. Then we have 

272 
The value of the atom-matter coupling parameter (37) essentially depends 
on the peculiarity of the atom-matter interaction. As the models of Section 4 
show, one should expect that a <^ 1. If the transition frequency of the doped 
atoms lies outside the polariton band gap, then the atom-matter resonance is 
possible, when Uks ~ ^o- Moving the atomic frequency into the gap makes such 
a resonance more and more difficult. Far inside the gap, where there are no 
elementary excitations of matter, this resonance becomes impossible. It follows 
from (53) that the relation between the atom-matter coupling parameter, aout, 
corresponding to the case when the atomic frequency is outside the gap, and the 
parameter ains, when the frequency is far inside the gap, is roughly speaking, as 

Oiout fi^ 

a- r2 ' 

for a sufficiently large polariton band gap. The decrease of a leads, as is clear 
from (80), to the increase of the delay time to. 

One can also notice that the coupling parameters g and a are not independent, 
but a depends on g. This dependence, for g ^ 1, is approximately as a ~ g^'^. 
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Therefore, strong atomic interactions diminish a, thus, increasing the delay time 
(80). During the interval < t < to, there is a temporal localization of light even 
for rather large parameters (?, such that g( > 1, but then the process of photon 
deconfinement starts. If to becomes comparable with Ti, one cannot omit in Eqs. 
(39), (40) the term containing 71. Then one should resort to numerical solution 
of these equations, which will be considered in a separate paper. 
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